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Abstract 

A finite difference formulation is presented 
for sound propagation in a two-dimensional straight 
soft-walled duct with uniform flow. The difference 
analysis is developed in terms of complex notation. 
The governing acoustic difference equations and the 
appropriate displacement boundary conditions asso- 
ciated with uniform flow are presented. Example 
calculations are presented for the sound attenuation 
in straight hard and soft-walled ducts. At present 
the finite Mach number case is solved only for the 
one-dimensional hard walled duct. 

I ntroduct i on 


In the aircraft industry, noise suppressors 
are now an integral part of turbofan engine instal- 
lation design. In particular, suppressors must be 
optimized for noise attenuation with minimum weight, 
length, and total pressure loss. Presently, esti- 
mates for the attenuation of acoustic energy (noise) 
in a suppressor are based on experimental data or on 
analytical studies' 1 "?) of acoustic waves propagat- 
ing into infinitely long straight ducts with uni- s 
form impedance (absorbers) along the walls. The 
analytical techniques, however, are limited to 
relatively simple geometries. 

There is a need for more flexible suppressor 
design techniques which can handle acoustical flow 
field complications arising from (1) axial varia- 
tions in Mach number, (2) axial variations in wall 
impedance associated with the effect of Mach number 
on the acoustical properties, or (3) sel f- genera ted 
noise sources at the boundary. To meet these needs, 
the present paper develops a numerical finite dif- 
ference technique for uniform mean flow. This tech- 
nique represents an extension of a difference theory 
developed for noise propagation in ducts without 
steady flow.' 8 ' The addition of a uniform Mach num- 
ber flow complicates the previous no-flow theory by 
the introduction of convective terms in both the 
governing equations and boundary conditions. 

The governing acoustic difference equations and 
the appropriate displacement boundary conditions as- 
sociated with uniform flow are presented. Immedi- 
ately following the mathematical development of the 
difference technique, numerical solutions are pre- 
sented for one and two dimensional hard and soft- 
walled ducts. The results are compared with the 
corresponding analytical results. 

List of Symbols 
c speed of sound 

dB decibels, Eq. (26) 

E acoustic power, Eq. (25) 

f frequency 

H channel height, see fig. 1 


I acoustic intensity 

i rrr 

L length of duct 

M Mach number, Jf/c 

m total number of grid rows 

n total number of grid columns 

p # pressure fluctuation, 

p’ (x 1 , y 1 , t) 

Pi Pi (*'. y 1 ) 

p dimensionless pressure fluctuation, 

P (x, y) = p # ’/p a 

P A amplitude of pressure fluctuation at duct 

entrance 

t time 

17 mean flow velocity i n X direction 

u dimensionless acoustical particle velocity 

in# direction 

AT dimensionless acoustic particle velocity 

y d i rect i on 

X dimensionless axial coordinate, x'/H 

axial grid spacing 

y dimensionless transverse coordinate, y'/H 

Ay transverse grid spacing 

Z acoustic impedance 

specific acoustic impedance, 
dimensionless frequency, Eq. (13) 

^ acoustic conductance ratio 

^ density 

0~ acoustic susceptance ratio 

&■ specific acoustic resistance 

^ specific acoustic reactance 

\jJ circular frequency 


Superscr i pts : 

J 1 or 2 

K 2 or 1 


J - 


1 



' prime indicates a dimensional quantity Equations (4) through (6) into Equation (1) through 

(3) and introduction of the following dimensionless 
* complex conjugate parameters 


(1) 

real part 

it 

4. 


(7) 

(2) 

imaginary part 

- * 


(8) 

Subscri 

e 

pts: 

exi t condi tion 

*- -f 


(9) 


i , j i axial index, j transverse index, 

see fig. 2 

w wall 

x axial position 

y transverse position 

Governing Equationa and Boundary Conditions 

The propagation of sound in a two dimensional 
rectangular duct, as shown in Figure la, is describ- 
ed by a solution of the continuity and momentum 
equations with the appropriate impedance boundary 
conditions. In this section, the basic governing 
equations and boundary conditions will first be non- 
d i mens i ona 1 i zed to show the important property 
groups. Next, the difference form of the governing 
equations and boundary conditions will be developed. 


yield the dimensionless steady state conservation 
equat i ons 


r-4fk# 

-7 


(10) 

Xi 

it 

p. 

t* 



(11) 


-igl® 


(12) 

di mensi on less 

frequency 

is given as 





(13) 


For practical linear design, )? is on the order of 
unity, which will lead to a manageable grid size in 
the numerical representation. 


Continuity and Momentum 


The linearized equations for mass and momentum 
conservation can be written for a cartesian coordi- 
nate system in the following form: 


(conti nu i ty) . , , „ 

(x-momentum) 


(y- momentum) 




-Hf 


(1) 

( 2 ) 
(3) 


The above dimensionless equations apply to the 
scaled cartesian coordinate system shown in Figure 
1(b). The dimensionless height y ranges between 
zero and 1 while the dimensionless length x ranges 
between zero and L/H. The solution of Equations 
(10) through (12) with the appropriate boundary con- 
ditions will lead to the determination of the duct 
attenuat i on. 

Wave Equation 

Equations (10) through (12) are now combined to 
yield the dimensionless wave equation for uniform 
flow 


The assumptions involved in the derivations of 
Equation (1) through (3) are given in most acoustic 
texts and will not be discussed herein. The single 
prime, 1 , is used to denote dimensional quantities. 
It will be removed when the equations are made di- 
mensionless. All the symbols used in this report 
are defined in the list of symbols. 


Equation (14) in difference form will be solved to 
determine the pressure in the duct. 

When using the exponential notation displayed 
in Equations (4) through (6), the dimensionless 
pressure and velocities have in general both real 
and imaginary parts. Thus, 


As customary for steady state, the solutions 
for the dependent pressures and velocities are as- 
sumed to be of the form 


(4) 


(5) 


(6) 

The final solution will be represented by 

the real 


part of the above quantities. Substitution of 


/Y*,y)= 

(15) 


( 16 ) 

s* (*<#)- 

(17) 

Consequently, Equation (14) can be broken into 
its real and imaginary parts by substituting Equa- 
tion (15) into Equation (14). The resulting two 
equations can be represented as a single equation 
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+ (rn> P ' D (18) 


A more general entrance condition where p(') de- 
pends on y is equally easy to treat. 

Exit Boundary Conditions 


where J equals 1 or 2 and K equals 2 or 1 , re- 
spectively. The two equations represented by Equa- 
tion (18) represent the basic governing equation 
for noise propagation which will be solved later by 
a difference formulation. As can be seen in Equa- 
tion (18), the real pressure p(') is coupled to the 
imaginary pressure p(2) through the term contain- 
ing Mach number to the first power. For zero Mach 
number, p(') and p(2) will be coupled by the im- 
pedance conditions along the walls. 

Whenever the J and K symbols are used lat- 
er in this paper, as in Equation (18), the single 
equation should be interpreted as representing two 
equations containing the real and imaginary contri- 
but i ons . 


Displacement Boundary Condition 

The boundary condition at the surface of a 
soft-wall duct requires that at the wall the acous- 
tic displacement of a particle in the fluid be equal 
to the displacement of a particle just inside the 
soft wall.v5) |n addition, the pressure and velocity 
fields at the boundary can now be related in terms 
of the specific acoustic impedance^ , where 

S AT 


As a consequence of these two relationships, it can 
be shewn that the appropriate boundary condition for 
slip flow is 



It is convenient to express the reciprocal of 
the impedance ratio J in terms of the conduc- 
tance ratio K and the acoustic susceptance ratio 
ff" , that is, 


Finally, the impedance must be specified at 
either the duct exit or at some point in the far 
field. The particular exit condition used in the 
example problem presented at the end of this paper 
will be discussed in a later section. 

Axial Acoustic Power 


The sound power which leaves a duct and reaches 
the far field is related to the axial intensity at 
the duct exit which for uniform flow is given in 
reference i as 


+ f c (u*"u. +/r** ir)J } 


(23) 


Substitution of Equations (15) through (17) into 
Equation (23), and non-d i mens i ona 1 i z i ng the result- 
i ng equat ion yei ids 


J-- Z&r I' 

fk 


]?v** '"W 

+j*.[r*>+ p'* 2 ] 


(24) 


*2 


The total dimensionless acoustic power is the integ- 
ral of the intensity across the test section 


By definition, the sound attenuation (the de- 
crease in decibels of the acoustic power from x = 0 
to x) can be written as 


-f- = X.- * r ( 2 °) 


where /j is the acoustic admittance ratio. 
Substituting Equations (15) and (20) into Equation 



Equation (21 represents the boundary condition on 
pressure to be used at the treated surfaces of the 
duct . 


AJB= loJpl-g] ( 26 ) 

Finite Difference Formulation 

instead of a continuous solution for pressure, 
the pressure will be determined at isolated grid 
points by means of the finite difference approxima- 
tions, as shown in Figure 2. This enables changing 
the differential equations to a system of algebraic 
equations for the pressure and velocity at each grid 
point. The governing equations and boundary condi- 
tions can be approximated in difference form(9) by 
using either a Taylor series expansion, a variation- 
al, or integral formulation. In this acoustic prob- 
lem, where the gradient is specified along a bound- 
ary, the integration method for generating the fi- 
nite difference approximations is most convenient. 


Entrance Conditions 


Wave Equation 


For the entrance pressure profile, the assump- 
tion used by Rice''' of a uniform profile will be 
used herein. At the present time, there is insuf- 
ficient information available to improve on this as- 
sumption. For a uniform pressure profile, as will 
be used in the example considered later, 




( 22 ) 


The wave equation in finite difference form is 
developed by applying the integration method (ref. 9, 
p. 168) to. the cell marked no. 1 in Figure 2. That 
is, the wave equation, equation (18), is integrated 
over the entire area of the cell marked number 1 in 
Figure 2. The details necessary for the development 
of the difference equations are given in reference 9 
as well as many other texts; consequently, only the 
final results are now presented. The difference 
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form of Equation (18) becomes: 


to* 


(3) 


(27) 


- (!-'#■) n«,j 

-(-ifzTTSlMAX 

+ (-/fzv/1 M AX P***. =. 0 


Displacement Boundary Condition 

Along the upper and lower boundaries, cell no. 
2, the wave equation, Equation (18), is again inte- 
grated, only this time Equation (21) is involved in 
the integration along the boundary. The difference 
equation that applies in cell no. 2 is 


to 

■a/ftl—t 


-t%K . 

<™K*£+£tr>,'rg 

+ -ZAfV 1 p' 7 * 

- [i (i-"*) - MKg (28) 

-r-ifrapm _* <r**+»'x t /‘ 

2 z¥)fAp *** 


M 


At I = -iztr/l p| 

**■* IJT^T r 'e 


(29) 


In terms of the real and imaginary pressures, this 
equation becomes 


P (J) ) = -(->y(zrrr>{) f (K} ) 

«)-% -4 d+M -4 


(K). 


(/+M 


(30) 


We wi 1 1 use Equation (30) as the boundary condition 
at the duct exit. Using Equation (30), the differ- 
ence representation of Equation (18) at the duct 
exit (cell no, 3) becomes 

2 I Aft V'ljJ 

+[>-** +<%f 
-4(3?) 

-r-lfziOlM AK 

+ £l) Z-rrjfAX P„ y j = D 


(31) 


-Hfteapai * C 

-0 

Ex i t Cond i t i on 

The exit condition to be presented now will al- 
low the numerical solutions to be compared to ana- 
lytical results for wave propagation in an infinite- 
ly long duct. The entrance region of length L in an 
infinitely long duct with uniform impedance will not 
have reflections at any position in the duct. The 
wave propagation in the entrance region of the infi- 
nite duct can be represented with a finite length of 
duct L by closing the exit impedance at L so 
that no reflection would occur. Consequently, if 
the numerical exit condition in the numerical analy- 
sis is chosen to eliminate reflections at the exit 
of the duct, the numerical and analytical results 
should be in close agreement. 

For a plane wave propagating in a hard-wall 
duct, the condition for no reflections at the duct 
exit is an exit impedance of fC which leads 
through the momentum equation to- the following exit 
cond i t i on 


Further details on the theory and limitations of 
Equation (31) are given in reference 8. 

Corner Condition 

At cell no. 4 in Figure 2, we apply the same 
integration technique and get a difference equation 
similar to the other equations presented. 

Acoustic Particle Velocity 

Recall from Equation (24) that the intensity 
is a function of the acoustic particle velocities. 
We must solve the two momentum equations, Equation 
(11) and (12) for the velocities. In terms of the 
real and imaginary terms, Equations (II) and (12) 
can be rewritten as 
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where the condition of i rrotat i ona I i ty 

iljgf s igg- 

was used in developing ^quati on (33)- 

Since the pressure gradients in the above equa- 
tions can be obtained directly from a numerical so- 
lution to the wave equation, Equations (32) and (33) 
can be treated as ordinary differential equations. 

As a result, a standard fourth-order Runge-Kutta 
integration is used for the solution of equation 
(32). The solution of Equation (33) can be found 
by a simple differentiation process' once u is 
known . 

The initial conditions for the Runge-Kutta so- 
lution can be found directly from an earlier part 
of the present work. In developing the displace- 
ment boundary condition, it also can be shown that 
the velocity at the wall is 

In a similar manner, the u velocity at the exit 
can be found directly from the exit impedance re- 
lations hi p 

Urn ~ Z7r *1 rfe (35) 


Axial Intensity 


In terms of the difference notation, the axial 
intensity as given by Equation (24) can be expressed 

as i r to to m to -| 

I hj 'Z-r* + Pc ’j Ui 'iJ 


*T] 

L -J 


+M 


2 (ZTT>1)* 
ft)* 


[“ 




■f U- ; 




tzi 
AT- • 


’] 


(36) 


The total power across a particular cross section, 
as given by Equation (25), is written in difference 
notation as 




4* -I 






(37) 


equations which can be expressed in matrix notation 
as 



where the A submatrix has a form typical of those 
matrices found in two-dimensional heat conduction 
problems, while the C matrix represents the coup- 
ling that occurs between the p' 1 ' and p(^) pressures 
through both the governing differential equations 
and impedance boundary conditions. The C matrix is 
a tridiagonal matrix. The F column vector contains 
the initial conditions. 

The frequency term (ZTTJtAK) 7 - in equa- 

tion ( 27 ) is a subtractive term in the main diagon- 
al element of the coefficient matrix; consequently, 
for sufficiently high frequency or spacing param- 
eters, the matrix will no longer be positive defi- 
nite (ref. 10, pg. 68). As a result, conventional 
iteration techniques can not be used. However, 
matrices of the form of Equation (38) can be solved 
by elimination techniques. In particular, the Gauss 
elimination technique will be used to find a solu- 
tion to the example problems which now follow. At 
present, the finite difference technique can not be 
applied to problems which require a large number of 
grid points because a subroutine based on a square 
full matrix is used for the solution. This subrou- 
tine requires large amounts of computer storage. 
Iteration techniques or more efficient closed form 
solutions will have to be developed to overcome 
present grid size limitations. 

Example Calculations 

Example calculations are presented in this sec- 
tion so that a comparison between the finite dif- 
ference solutions and corresponding analytical solu- 
tions can be made. First, some solutions for wave 
propagation with no mean flow (zero Mach number) are 
presented. Next, solutions for wave propagation in 
a hard wall duct with uniform flow are presented. 

Zero Mach Number 


Hard wa 1 1 duct . Numerical and analytical val- 
ues of the pressure are computed for the case of a 
one dimensional hard wall duct with a fC exit im- 
pedance. Because there is no variation of pressure 
in the y direction, the two dimensional grid lattice 
shown in figure 2 reduces to a one- di mens i ona 1 lat- 
tice as shown in the upper sketch of figure 3- The 
calculation was made for a hard wall (Z = o© ) wi th 
an L/H of lj an inlet plane wave with a dimension- 
less frequency^ equal to I, and zero Mach number. 
The analytical and numerical values of the acoustic 
pressure profiles along the duct are shown in Fig- 
ure 3- As seen in Fiqure 3, agreement between the 
numerical and analytical results is good. 


By evaluating &C at the entrance and exit 
positions, taking the log of their ratio and multi- 
plying by 10, as indicated in Equation (26), the 
sound attenuation for the duct is determined. 

Matrix Solution 

The collection of the various difference equa- 
tions at each grid point forms a set of simultaneous 


Soft wall duct ■ As another example of the 
finite difference formulation, the noise attenuation 
will now be calculated for a two-dimensional duct 
with a dimensionless frequency 7^ of 0.6, an L/H of 
0.5, and a finite value of wall impedance. At the 
present time for two dimensional solutions, the 
technique is limited to L/H values of 05 or less 
because of matrix size limitations discussed pre- 
v i ous 1 y . 
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Sound attenuation values were obtained by cal- 
culating the appropriate values of and 0 ~ from 
Equation (20), solving the matrix equation ( 38 ) by 
Gauss elimination for the pressure distribution, 
solving Equations (32) and (33) for the velocity 
distribution, and finally calculating the sound 
attenuation at the exit from Equations (37) and (26) . 

Because of the matrix size limitations in the 
subprogram used to solve Equation ( 38 ), the tech- 
nique of varying grid size was used. This is 

illustrated in figure 4 where sound power attenua- 
tion is plotted against grid spacing 4 y. The true 
attenuation is found by extrapolating the grid size 
to zero (dotted portion of the curve). As seen in 
Figure 4, the extrapolated numerical values for the 
zero Mach number case are in good agreement with the 
analytical values calculated from theory similar to 
that presented by Rice in reference 1. 

The calculations shown in Figure 4 are based 
upon a four point numerical derivative formula rath- 
er than the two point formula used in reference 8. 

As a consequence, the grid point extrapolation tech- 
nique for equal 0.6 is linear. The linear ex- 
trapolation more clearly defines the answer for zero 
grid size, than the two point formula used in refer- 
ence 8. 

Finite Mach Number 

Hard wal I duct . The numerical and analytical 
values of both the pressure and velocity are com- 
puted for the case of a one dimensional hard wall 
duct as was considered before (Figure 3), but with 
Mach numbers of 0-5 and -0.5 . The results are shown 
in Figures 5 and 6. As seen in both of these fig- 
ures, the agreement between the numerical and ana- 
lytical results is good. In Figure 5, the positive 
Mach number lengthens the effective wavelength of 
the sound, while in figure 6, the negative Mach 
number compresses the effective wavelength compared 
to the no-flow case (Figure 3)- 

Soft wal 1 duct . Currently, the difference 
equations for the two dimensional soft walled duct 
with finite Mach number are being programmed and 
debugged. At the present time, some numerical 
stability problems are occurring at low values of© - 
and values of L/H near 1. We are investigating the 
possibility that a finer grid structure will elimi- 
nate these instabilities. 

Concluding Remarks 

A finite difference theory for sound propaga- 
tion in a two-dimensional, soft-walled duct with 
uniform mean flow has been presented. The governing 
equations and boundary conditions are presented in 
difference form for both the real and imaginary 
pressure and velocity. 

The numerical theory is shown to be in good 
agreement with the corresponding exact analytical 
results. However, because the solution matrix for 
the acoustical flow field is not posi tive def ini te, 
conventional iteration techniques cannot be used to 
solve the difference equations. Before the finite 
difference formulation can be applied to more com- 
plicated problems which require a large number of 
grid points, iteration techniques or more efficient 
closed form solutions will have to be developed to 
overcome present grid size limitations. 


The finite difference formulation is flexible 
and should be a powerful tool in studies of inlet 
and exhaust ducts of turbofan engines. 

References 

1. Rice, E. J.: Performance of Noise Suppressors 

for a Full-Scale Fan for Turbofan . Engi nes , NASA 
TM X-52941 , Jan. 1971 ■ 

2. Eversman, W. : The Effect of Mach Number on the 

Tuning of an Acoustic Lining in a Flow Duct. 

J. Acoust. Soc. Amer. 48, No. 2 (Part l), pp. 
425-428, 1970. 

3- Lambert, R. F.; and Tack, D. H.: Influence of 

Shear Flow on Sound Attenuation in a Lined Duct. 
J. Acoust. Soc. Amer. 38 , 655-666 (1965). 

4. Pr i dmore- Brown, D. C.: Sound Propagation in a 

Fluid Flowing Through an Attenuation Duct. J. 
Fluid Mechanics 4, 393-406 (1958). 

5- Mungur, P. ; and Gladwell, G. M. L.: Acoustic 

Wave Propagation in a Sheared Fluid Contained 
in a, Duct. J. Sound Vib. 9 (I), 28-48 (1969)- 

6. Savkar, S. D-: Propagation of Sound in Ducts 

with Shear Flow. Journal of Sound and Vibration 
19, 355-372 (1971). 

7- Shankar, P. M.: On Acoustic Refraction by Duct 

Shear Layers. Journal of Fluid Mechanics 47, 
81-91 (1971)- 

8. Baumeister, Kenneth J.: Application of Finite 

Difference Techniques to Noise Propagation in 
Jet Engine Ducts, NASA TM X-68261 , Paper to be 
presented at ASME Winter Annual Meeting, Detroit, 
Nov. 1 973 - 

9. Varga, Richard S.: Matrix Iteration Analyses, 

Prentice-Hall, Inc., 1962. 

10. Morfey, C. L-: Rotating Pressure Patterns in 

Ducts: Their Generation and Transmission. J. 

Sound Vib. (1964), I, 60-87- 

11. Morfey, C. L.: A Note on the Radiation Effic- 

iency of Acoustic Duct Modes, J. Sound Vib. 

(1969) 9(3), 367-372. 

12. Moon, Parry; and Spencer, Domina Eberle: Field 

Theory for Engineers. Van Nostrand, 1 96 1 . 


6 



-7644 



(A) DIMENSIONAL CARTESIAN COORDINATE (B) DIMENSIONLESS (SCALED.) 
SYSTEM, p', v', u\ u). CARTESIAN COORDINATE 

SYSTEM, p, u, v, 77 . 

Figure 1. - Schematic of two dimensional suppressor duct. 



Figure 2. - Coordinate and grid point representation of two dimensional 
soft walled duct. 
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Figure 3. - Analytical and numerical pressure profiles for one 
dimensional sound propagation in hard wall duct for dimen- 
sionless frequency 77 = 1 and M =0. 
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AV ANALYTICAL SOLUTION 
RICE (REF 1) 



7j=0.6 (DIMENSIONLESS FREQ) 
n = 5 (AXIAL GRID POINTS) 
0=0.160 (SPECIFIC RESISTANCE) 
X = -0. 34 (SPECIFIC REACTANCE) 
WHERE £ = 0 + ix 


/-EXTRAPOLATED 



77 = 1 
n = 10 
8 = 0. 230 


X-0.55 


(DIMENSIONLESS FREQ) 
(AXIAL GRID POINTS) 
(SPECIFIC RESISTANCE) 
(SPECIFIC REACTANCE) 


2L 

0 



.3 


Ay 


Figure 4. - Effect of spacing on attenuation in two dimen- 
sional duct for L/H=0. 5 and M=0. 
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Figure 5. - Analytical and numerical 
pressure and velocity profiles for 
one dimensional sound propagation 
in hard wall duct for dimensionless 
frequency 17 = 1 and M =0. 5. 
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Figure 6. - Analytical and numerical 
pressure and velocity profiles for 
one dimensional sound propagation 
in hard wall ducts for dimension- 
less frequency 77 = 1 and 
M = -0. 5. 
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